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I use the two-step density-matrix renormalization group method based on two-leg ladder expan- 
sion to show numerical evidence of a plaquette ground state for J2 — 1.3Ji in the Shastry-Sutherland 
model. I argue that the DMRG method is very efficient in the strong frustration regime of two- 
dimensional spin models where a spin-Peierls ground state is expected to occur. It is thus comple- 
mentary to quantum Monte Carlo algorithms, which are known to work well in the small frustration 
regime but which are plagued by the sign problem in the strong frustration regime. 



A number of studies have been devoted to the Shastry- 
Sutherland model (SSM) [H, i, i, i, i, S, 0, i, i, E [ll 
12|. This interest is motivated by the relevance of the 



SSM to the physics of the two-dimensional spin gap sys- 
tem 5'rC'U2 (503)2 HI. The SSM is a frustrated antifer- 
romagnetic model on a square lattice whose Hamiltonian 
is written as 



H — Ji SjSj -I- J2 SiSj, 



(1) 



<«j> 



[i,3] 



where < i,j > represents nearest-neighbor sites and 
stands for the next-nearest neighbors on every other pla- 
quette in the pattern shown in Fig[T] There is a general 
agreement that in the weak frustration regime Ji 3> J2, 
the SSM is Neel ordered while in the strong frustration 
regime Ji <C J2, the model is a valence bond solid. In 
fact, Shastry and Sutherland showed that for J2 > 2Ji 
the wavefunction made of the product of local orthog- 
onal dimcrs is an exact eigenstate. Numerical simula- 
tions based on series expansions and exact diagonaliza- 
tion have pushed the dimer phase boundary down to 
J2 ~ I.5J1. The estimated boundary of the Neel phase 
is J2 « 1.2Ji. 

However, the nature of the ground state for 1.2Ji < 
J2 ^ I.5J1 so far has remained controversial. A mean- 
field Schwinger boson approach found that the interme- 
diate phase is helical 0]. Perturbation theory Q and 
series expansion Q studies predicted a direct first or- 
der transition between the Neel and the dimer phases. 
A large N study [rj also predicted a helical phase. In 
addition it suggested that a broader phase diagram con- 
tains a plaquette phase which might occur if fluctuactions 
were included. A subsequent series expansion analysis Q 
predicted a plaquette phase with a spin gap. This con- 
clusion was criticized in another series expansion study 
which suggested a possible gapless phase whose nature 
was unclear Q. The existence of the intermediate phase 
was also suggested in a renormalization group analysis 
An exact diagonalization study on a = 32 site 
system has concluded to a plaquette phase. The Monte 
Carlo method, which is very effective for spin systems in 
absence of frustration , is plagued by the sign problem 
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FIG. 1: (a) Bond patterns in the Shastry-Sutherland lattice. 
The Ji bonds are between the nearest neighbors. The J2 
bonds are along every other diagonal, (b) the two-leg lad- 
der: starting point of the two-step DMRG approach to the 
two-dimensional lattice. (c) Orthogonal dimers: starting point 
more adapted to the dimerized phase. It is obtained by cov- 
ering the two-dimensional lattice allong the diagonals. 



in this regime of strong frustration. Knowledge of the ex- 
act ground state phases is essential; this could serve as a 
starting point in variational investigations of the nature 
of superconductivity that might arise upon doping [T^. 

In this letter I present numerical evidence of the pla- 
quette phase at J2 = 1.3Ji. For this purpose, I will 
use the two-step density- matrix renormalization group 
(DMRG) method [11 [ll]. The DMRG [H has provided 
a breakthrough in the study of quantum Hamiltonians 
in one dimension. Extensions of the DMRG to two- 
dimensional Hamiltonians have been less effective. Liang 
and Pang 19] found that as the linear dimensions of the 
system grow, the number of the reduced density matrix 
states needed to maintain accuracy grows exponentially. 
This problem is particularly severe for quantum antifer- 
romagnets in their ordered phase. The spontaneous sym- 
metry breaking which takes place in the thermodynamic 
limit is due to the collapse of an infinite number of ex- 
cited states onto the ground state. The implication for 
finite systems, in the parameter regime where long-range 
order occurs, is a near degeneracy of a large number of 
low-lying states with the ground state. Each state within 
this large set would carry the same weight in the reduced 
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FIG. 2: Gaps as a function of L: for a single ladder for J2 = 
(circles), J2 = 1.3Ji (squares), and for two-dimensional 
systems for J2 = (diamonds), J2 = 1.3Ji (triangles). 



FIG. 3: Edge-to-center correlation as a function of Ly-. for a 
single ladder for J2 =0 (circles), J2 = 1.3Ji (squares), and for 
two-dimensional systems for J2 = Q (diamonds), J2 = 1.3Ji 
(triangles) . 



density matrix. For this reason, standard DMRG sim- 
ulations of spin Hamiltonians are limited to systems of 
about ten sites wide. 

Recent developments by the author 14 , [3, Ell have 
shown that the DMRG could be very useful for 2D models 
in the region of the parameter space where this techni- 
cal difficulty is less severe or even absent. This occurs 
for instance for spatially anisotropic antiferromagnets or 
in the highly frustrated regime of isotropic magnets. In 
this latter case, general arguments from the large N ap- 
proach [l3| suggest that the ground state is a collection of 
weakly-coupled dimers or plaquettes. Presumably deep 
in the disordered phase, because of the presence of a spin 
gap, the ground state would be dominant in the reduced 
density matrix. This is more favorable to a DMRG sim- 
ulation. It is usually in this regime that the sign problem 
is most severe in QMC simulations. Hence the DMRG 
would be complementary to QMC for frustrated models 
with disordered phases with broken translational sym- 
metry. The excited states would become more and more 
important as the coupling is moved toward the bound- 
ary with the magnetically ordered phase. An approach 
based of these ideas recently has been applied to find the 
gro und-state phase diagram of the checkerboard model 
[1^. The same approach is applied here to the SSM. 

I start with a single two-leg ladder with L rungs as 
shown in Fig[TJ By doing so, I implicitly assume that 
the inter-ladder interactions are small. Strictly, this is 
true only in the magnetically disordered phase and will 
be justified a posteriori. However, such a starting point 
could also be justified qualitatively even for the magnet- 
ically ordered phase where inter-ladder interactions are 
not small. In the valence bond represensation of quan- 
tum antifferomagnets, the wave function is written as 



a (y)e{a} 



where |(ij)) = (| tDij - | iT)ii)/V2 is a dimer wave func- 
tion between the sites i and j, and {a} is a configuration 
of dimers. Such valence bond representations are quali- 
tatively good for both the disordered and ordered phases 
[2^ . In the present approach, two- leg ladders are the 
building blocks. The wave function is written as 



^ ^ ^ladders 
ladders 



n 

ladders 



ladder 1 



(3) 



(2) 



where ^ladder is an eigenfunction of a single ladder 
Hamiltonian. Given that the lowest ^ladder is dominated 
by a product of dimers, 4" bears some similarity with ^. 
However, the set of ^ladder ^ includes excited states on 
the ladder, the structure of 4" is thus much more com- 
plex than a simple short-range dimer expansion. When 
the ground state is made of weakly-coupled plaquettes 
or dimers, it would be expected that this representation 
would yield quantitatively good results as well. But the 
essential point is that this approximation does not neces- 
sarily assume that the ground state is disordered. It will 
be shown below that a magnetically ordered state can be 
reached as well, though with less accuracy than in the 
disordered phase. 

The results for an isolated two-leg ladder were obtained 
for J2 = and for J2 = 1.3Ji. The conventional DMRG 
is known to yield highly accurate results for the ground 
state and the low- lying states 2l|. A set of the low- 
lying ^ladder is obtained by targeting the spin sectors 
from 5* = to 5 = ±4 and by keeping up to ra = 144 
states. This is enough to maintain the truncation er- 
ror below 10~^ in all cases. There is a spin gap in the 
thermodynamic limit as in two cases seen in Figl2l The 
center-to-end correlation functions Cl = (8^/2+1, rSi.r), 
where r = 1,2 is the leg index, shown in Fig[3] decay 
exponentially in both cases. But the short-range corre- 
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FIG. 4: (a) Ground state energy of two-leg ladders: for the 
Shastry-Sutherland model at J2 = 1.3Ji (circles) and for the 
checkerboard model at J2 = Ji (squares), (b) Short-range 
correlations for the Shastry-Sutherland two-leg ladder with 
Lx ~ 24: for J2 = 0, Ci^ ^ (circles), Cr (squares); for J2 = 
I.3J1, Ci^ ^ (diamonds), Cr (triangles). 



lations in Fig|4] reveal that J2 = and J2 = 1.3Ji be- 
long to two different phases of the ladder. When J2 — 0, 
the dominant short-range correlation are the rung dimers 
Cr{i) — (Si,iSi.2). Cr is stronger than the correlations 
along a leg C;^ ^(i) = {Si,rSi+i^r)n,d for plaquettes with 
no diagonal bond (n) or with a diagonal bond (d) ; when 
J2 = 0, — Ci^. Both Cr{i) and Ci^^{i) are inde- 
pendent of i, except for small variations at the bound- 
ary. For J2 = 1.3Ji, |Ci„| > \CiJ and the bond pattern 
shows strong alternations as function of i as seen in FigU) 
At the same time, Cr ~ Ci^ indicating that the system 
is now in the plaquette phase. The plaquette-plaquette 
interaction is given by |QJ. It is not negligible as in 
the checkerboard ladder It is about one third of 

the intra-plaquette interaction in the Shastry-Sutherland 
ladder, while it is only 5% of the intra-plaquette in the 
checkerboard ladder at the isotropic point J2 = Ji . This 
is also seen in the ground state energy Eg of the lad- 
ders shown in FigH) In the checkerboard ladder, there 
is a relatively small renormalization of Eg = — 0.5086 Ji, 
which is not very far from Eg — —0.5000 Ji of an isolated 
plaquette. The renormalization is more important in the 
Shastry-Sutherland ladder where Eg — — 0.5263 Ji. Nev- 
ertheless, the value of Ci^ shows that even the SSM is in 
the weak-coupling regime of plaquette-plaquette interac- 
tion. 

The second step of the two-step DMRG consists in 
projecting the Hamiltonian ([1]) onto the basis states of 
the tensor product of ^ladder s and solving the result- 
ing effective Hamiltonian, which is one-dimensional ( in 
the transverse direction) with the usual DMRG. Strictly 
speaking, this can only rigorously be justified if the inter- 
ladder coupling is small. Apparently this is not the case 
for Hamiltonian ([T]), since neither Ji nor J2 is always 
small. But as seen above, effective small interactions can 
be generated by the inherent competition between Ji and 
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FIG. 5: Irrelevant (squares, J2 = 1.3Ji) and relevant (circles, 
J2 = 0) flows from the single ladder as function of Ly : (a) for 
the spin gap, (b) for the center-to-edge correlation function. 



J2. The coupling between ladders in the SSM is also given 
by |Ci_j|, as can be seen in FigH] In Ref [iBt .it has been 
shown that when the ratio, p = SE/J^ff between the 
bandwidth of the states kept, and the effective transverse 
couphng, which is Jeff « \CiJJi, is large {p > 4), the 
two-step DMRG yields results which are comparable to 
those of QMC. Typical values of p for the largest systems 
studied are 5 for J2 = 1.3Ji, and 2 for J2 = 0. This jus- 
tifies the isolated ladder starting point for J2 = 1.3Ji. I 
study lattices with L^xLy — Lx{L+l) = 8x9 to 24x25. 
Analysis of the performance of this app roach have been 
discussed in previous pubHcations [ij When the 

inter-ladder is turned on, the flows as the function of the 
number of ladders of A and Cl for J2 = and J2 = 1.3Ji, 
shown in Figl5]for L x (L-l-1) = 24 x 25, are very different. 

For J2 = 0, as expected from the existence of long- 
range order, A decays rapidly as Ly increases. A = in 
the thermodynamic limit as seen in Figl2] At the same 
time, Cl grows away from its small value found in the 
two-leg ladder. Cl becomes finite in the thermodynamic 
limit. The extrapolated order parameter M = \fC^ is 
found to be M = 0.1738. This is somewhat lower than 
the QMC results [111 M = 0.3070. Part of this discrep- 
ancy is due to the use of open boundary conditions which 
yield an undervaluated Cl- Better extrapolations can be 
obtained if the lattice sizes are reduced, and if periodic 
boundary condition are used. This relatively poor per- 
formance of the DMRG deep in the magnetically ordered 
phase is a consequence of the fundamental limitations of 
the DMRG when faced with an exponentially dense low 
energy spectrum, as discussed in the introduction. Nev- 
ertheless, this results shows that the approximation by 
of the exact wave function retains the correct qualitative 
behavior. Hence the DMRG could still be very useful in 
the magnetic regime as well. 

However, when J2 = 1.3Ji, the inter-ladder coupling 
does not qualitatively affect the physics of a single lad- 
der which is itself that of nearly isolated plaquettes as 
seen above. The irrelevant flows for A and Cl with 
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the number of ladders are shown in FigO A for the 
two-dimensional system is renormalized by about 20% 
from its single ladder value. The extrapolated gap A = 
0.4758 Ji is lower than A — 0.67Ji found in the checker- 
board model at the isotropic point J^. This is consistent 
with the fact that inter-plaquette interactions are more 
important in the SSM. For large lattices, Cl for the two- 
dimensional systems is practically identical to its ladder 
value. This suggests that the correlation length is very 
short. In this regime short-ranged plaquettes are domi- 
nant in the exact wave-function, hence 5* is an excellent 
variational approximation. 

Recently a variational wave function based on doping 
an orthogonal dimer wave function has been used to ex- 
plore the nature of an eventual superconductive state in 
SrCu2{BOz)2 [13]. The estimated value of the couplings 
in this compound are Ji « 85i^ and J2 = bAK Q . This 
places it at the boundary of the orthogonal dimer phase, 
not very far for the plaquette phase. Given the possible 
incertitude in this estimation and the fact that these val- 
ues may be affected by doping, it is worth exploring su- 
perconductivity upon doping the plaquette ground state 
as well. 

To summarize, I have argued that the DMRG tech- 
nique is a natural approach to study spin-Peierls phases 
that spontaneously arise in frustrated quantum antifer- 
romagnets. I have used the two-step DMRG to confirm 
the nature of the controversial phase between the Neel 
and dimer phases. The phase diagram of the SSM bears 
some similarity to that of the checkerboard model . In 
the checkerboard model, an additional Neel phase was re- 
cently found between the plaquette and the crossed dimer 



phases [16|, l25| . It is quite possible that this additional 
phase exists in the SSM in the vicinity of J2 ~ 1.5Ji. Un- 
fortunately, in the SSM, unlike the checkerboard model, 
the starting ladder does not have the full symmetry of 
the bond pattern of the 2D lattice. As the dimer phase 
is approached, the variational wave function used in this 
study is not optimal, as it leaves half of the spins un- 
paired in the orthogonal dimer phase. It will be more 
advantageous to start with the orthogonal dimer pattern 
shown in Fig[TJ 
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